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Abstract. We investigate the quasiparticle band structure of anatase Ti02, a wide 
gap semiconductor widely employed in photovoltaics and photocatalysis. We obtain 
GW quasiparticle energies starting from density-functional theory (DFT) calculations 
including Hubbard U corrections. Using a simple iterative procedure we determine the 
value of the Hubbard parameter yielding a vanishing quasiparticle correction to the 
fundamental band gap of anatase Ti02. The band gap (3.3 eV) calculated using this 
optimal Hubbard parameter is smaller than the value obtained by applying many-body 
perturbation theory to standard DFT eigenstates and eigenvalues (3.7 eV). We extend 
our analysis to the rutile polymorph of Ti02 and reach similar conclusions. Our work 
highlights the role of the starting non-interacting Hamiltonian in the calculation of 
GW quasiparticle energies in Ti02, and suggests an optimal Hubbard parameter for 
future calculations. 
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The wide gap semiconductor Ti02 is the subject of a vast research effort owing to its 
importance to many areas of technology, ranging from photovoltaics and photocatalysis 
to sensing devices, protective coatings and batteries p. As an example, nanostructured 
anatase Ti02 films sensitized with molecular dyes are the fundamental building blocks 
of dye-sensitized solar cells [2]. In order to understand and optimize the photovoltaic 
and photocatalytic properties of Ti02 it is important to complement experimental 
spectroscopies with reliable calculations of electronic excitations in this material and 
its interfaces. 

Band structures calculated within standard DFT underestimate the fundamental 
band gap of anatase Ti02 by about 1 eV. The underestimation of the band gap is 
accompanied by errors of similar magnitude in the individual addition and removal 
energies, and raises concerns about the calculation of energy-level alignments and charge 
injection dynamics at interfaces involving Ti02 [3]. 

In this context accurate quasiparticle methods for studying electronic excitations, 
such as the GW approximation in its various implementations |11 El El [3 IH] , still face two 
limitations, (i) On the one hand these methods are computationally demanding, and 
cannot be used yet for complex systems such as defects and interfaces. For example the 
minimal atomistic models required to describe donor /acceptor interfaces for solar cells 
involve several hundreds of atoms [H [TOl [11] , which are beyond the reach of current GW 
implementations, (ii) On the other hand, there is the issue of whether GW calculations 
should be performed at the GqWq level, or include some form of self-consistency [H [12] 
or off-diagonal matrix elements of the self-energy [12]. 

In this work we attempt to address the latter point by combining DFT+[/ and GqWq 
calculations, in the spirit of Kioupakis et al [13]. The rationale for this approach is that 
the use of DFT+[/ as the noninteracting Hamiltonian should provide a better starting 
point for perturbation theory, thereby legitimizing the use of the GqWq approximation. 
The Hubbard parameter is determined in such a way that a GqWq calculation starting 
from the DFT+^7 band structure yields a vanishing correction to the fundamental band 
gap. This is achieved through an iterative procedure consisting of a sequence of GqWq 
and DFT+[/ calculations, and can be seen as a simple and inexpensive way to mimic 
self-consistent GW . Using this procedure we obtain a band gap of 3.3 eV for anatase 
Ti02. When comparing the band structures obtained from DFT+t/ with or without 
quasiparticle corrections, we find that they are in reasonable agreement only near the 
band edges, and that the individual shifts of the bands determined from DFT+[/ are 
always larger than the corresponding quasiparticle shifts. Interestingly we find that the 
optimal value of U thus determined is transferable to the rutile polymorph of Ti02. 

We perform DFT calculations using the generalized gradient approximation of 
Perdew et al [Hj and norm- conserving pseudopotentials [15]. We explicitly include 
the 3s and 3p semicore states of Ti in our calculations. The electronic wavefunctions 
are expanded in a plane waves basis set with a kinetic energy cutoff of 200 Ry, and 
the Brillouin zone of Ti02 is sampled using a 6x6x6 Monkhorst-Pack mesh. All DFT 
calculations are performed using the Quantum ESPRESSO software package [16]. For 
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Figure 1. Convergence of GWU calculations with the energy denominator cutoff in 
the sum-over-states expansions of the polarizability and the Green's function of anatase 
Ti02: X2Av-^2bc gap (circles, left scale), and quasiparticle correction to the ^241, state 
(squares, right scale). The quasiparticle corrections were obtained by starting from a 
DFT+f/ calculation with U—7.5 eV. The lines are guides to the eye. 

definiteness we set the lattice parameters and atomic positions within the unit cell to 
the experimental values [l7j. In order to perform DFT+U calculations |18j we adopt 
the simplified rotationally invariant formulation of Cococcioni and de Gironcoli [19] , with 
the Hubbard-like corrections included using the atomic Ti 3d pseudo-wavefunctions for 
the projectors. It is worth stressing that, as in every DFT+[/ calculation, the value 
of U determined in this work is linked to the particular choice of the projectors [20] . 
Quasiparticle GW calculations are performed using the scheme of Hybertsen and Louie 
[5j. The frequency dependence of the screened Coulomb interaction is described using 
the Godby-Needs plasmon-pole model [21] , by evaluating the static polarizability as well 
as the polarizability at the imaginary frequency of 23 eV. For the screened Coulomb 
interaction we use a uniform and unshifted 4x4x4 Brillouin zone mesh. The exchange 
and the correlation components of the self-energy are described using plane waves basis 
sets with kinetic energy cutoffs of 80 Ry and 10 Ry, respectively. In the expansion 
over conduction states we include 568 unoccupied bands for calculating both the 
polarizability and the Green's function. The energy denominator cutoff, i.e. the energy 
of the topmost band used in this expansion referred to the valence band maximum, is 
16.6 Ry. Figure [U shows the convergence of the fundamental gap with the denominator 
cutoff. Increasing the denominator cutoff from 10.6 Ry (296 conduction bands) to 
16.6 Ry (568 bands) reduces the band gap by 0.02 eV. The individual corrections to the 
band edges converge more slowly, as shown in Figure [1] for the highest occupied state 
at the X point. All the GW calculations are performed using the SaX code [22j. 

The calculation of GW quasiparticle energies using DFT+^7 as the starting non- 
interacting Hamiltonian is performed as follows. We calculate the GqWq self-energy E 
using eigenstates and eigenvalues obtained from a DFT+^7 calculation. Then, for each 
Kohn-Sham state -0 of energy e, the quasiparticle energy E is obtained from: 

E = e + Z(e) (V'l S(e) - Kc - Vu 1^) , (1) 

J Our optimized lattice parameters are a=3.83 A and c—9.84 A and the internal parameter is u=0.206 
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Figure 2. Band structure of anatase Ti02 calculated using various methods. The 
gray lines represent standard DFT calculations [13], the black lines are obtained from 
DFT+[/ calculations [19] using our optimal Hubbard parameter J7=7.5 eV, and the 
blue disks are from GW calculations using DFT+C/ as the non-interacting Hamiltonian. 
In each case the top of the valence band is arbitrarily aligned with the zero of the energy 
axis. 

where Z is the quasiparticle renormahzation, 14c is the DFT exchange and correlation 
potential, and Vu is the Hubbard term in the Kohn-Sham Hamiltonian. The Hubbard 
contribution to the quasiparticle shift is conveniently evaluated as: 

(^1 Vu \^) = t- (^1 ^0 1^) , (2) 

where Hq is the Kohn-Sham Hamiltonian without the Hubbard term (but with the 
DFT+t/ density). In the following we will refer to this procedure simply as GWU 
method. We have compared the present approach with the method described by Jiang 
et al [23], where the dynamical contributions are evaluated at the energy e — (■?/'| V(/ 
and found the differences to be negligible. In the case of [/=0 equation ([T]) yields the 
standard GW quasiparticle energies obtained by using perturbation theory on top of 
DFT. 

The use of DFT+f/ as the starting Hamiltonian for the perturbation expansion 
introduces an additional arbitrariness in the procedure, relating to the choice of the 
Hubbard parameter U. Since we are using GW within the framework of perturbation 
theory, the natural choice for the starting Hamiltonian (and hence the value of U) is the 
one leading to the smallest quasiparticle corrections. This observation can be formalized 
into a criterion for choosing the Hubbard parameter: our optimal U is determined by the 
condition that, upon performing a GWU calculation, the fundamental band gap does 
not change. This approach has been demonstrated for the case of bcc hydrogen [13]. The 
optimal Hubbard parameter is determined by using an iterative procedure consisting of 
successive GWU and DFT+t/ calculations, until self-consistency is achieved. 

Figure [2] shows that the band gap of anatase Ti02 is indirect, with the top of the 
valence band located near the X point of the Brillouin zone, and the bottom of the 
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conduction band at the F point. Since the top of the valence band hes very close in 
energy with the states at the X point (<60 meV), for convenience in our analysis we 
take the states at the X point as representative of the valence band top. We denote 
the band extrema as ^24^ and r25c, respectively [21]. The values of the X2iv-^25c gap 
calculated within DFT+[/ for a range of values of U between and 10 eV are shown 
in Fig. [3l and are in good agreement with previous work [25] . For U =0 we recover the 
standard DFT band gap of 2.08 eV (indicated by A in figure [3]). 

Now we go through the iterative calculation of the band gap. We calculate the GW 
quasiparticle corrections starting from the DFT eigenstates and eigenvalues (i.e. from 
point A in figure |3]), and obtain a corrected band gap of 3.70 eV (point B in figure E]). 
The latter value is in good agreement with previously published works using the same 
methodology, 3.6-3.8 eV [Ml ESI EZl [28]. We note that the Godby- Needs plasmon- 
pole model [21j used in our calculations produces a band gap in good agreement (only 
0.08 eV larger) with calculations based on the contour deformation method [26], and 
does not suffer from the overestimation of the band gap reported for other plasmon- 
pole models [26]. The success of the Godby-Needs approach in reproducing contour 
deformation results has recently also been demonstrated for ZnO [30]. From figure [3] we 
see that the same band gap as in B can be obtained by using DFT+[/ and the Hubbard 
parameter [7=9.5 eV (point G). By repeating the previous procedure and applying GW 
corrections to the DFT+[/ eigenstates and eigenvalues obtained with [7=9.5 eV, we find 
that the quasiparticle corrections slightly reduce the band gap to 3.20 eV (point D in 
figure |3]). We iterate this procedure until the quasiparticle correction to the band gap 
becomes negligible. The iterations generate the converging series of Hubbard parameters 
0.0^9.5^7.0^7.5 eV, and the series of band gaps 2.08^3.70^3.20^3.29^3.27 eV. 
By the fourth iteration the DFT+[7 and GWU band gaps are both 3.27 eV and the 
Hubbard parameter is U=7.5 eV. This value of the Hubbard parameter is compatible 
with the range of values 2-8 eV reported in the literature for oxides containing Ti, 
cf. Nolan et al [29] and references therein. It is interesting to check whether the optimal 
Hubbard parameter [7opt=7.5 eV here determined for anatase Ti02 is transferable, i.e. 
(i) whether it generates a band structure in agreement with GWU, and (ii) whether it 
can also be used for other Ti02 polymorphs. 

In order to address the first point we compare in Fig. [2] the band structures 
calculated using DFT+[7 and those calculated using the GWU method, both with 
U=Uopt- Near the band extrema the two methods yield similar results, although 
differences in energy of up to 0.3 eV are observed for several bands. We also estimate that 
the GWU calculation yields a valence band width which is 0.3 eV larger than the starting 
DFT+[7 band structure. These differences are rather small when compared with the 
broadening of the valence photoemission spectra reported by Li et al [31] , corresponding 
to ~1 eV. While the band gaps and band width obtained using DFT+U and GWU are 
in reasonable agreement, tabled] shows that the corrections to the individual eigenvalues 
differ significantly. In particular, DFT+[7 consistently yields energies which are ~0.6 eV 
higher than in GWU. 
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Figure 3. Fundamental band gap (^24t>-r25c) of anatase Ti02 calculated using 
DFT+i7 (black squares), and GWU (blue circles), respectively. The lines are guides 
to the eye. 

In order to address the second point we consider the rutile structure. In our DFT 
and GWU calculations the fundamental gap is direct at the T point, therefore in tabled] 
we only report our results for the T24,v and T25C states. We note, however, that the 
energies of the conduction band minima of rutile at F, M, and R are within 60 meV 
from each other, and that the location of the minimum is dependent upon the calculation 
method. Given the small energy difference between direct and indirect gaps in rutile, in 
order to establish the precise nature of the gap it would be necessary to investigate the 
magnitude of the phonon- induced renormalization [32l [33l [341 [35]. Table [T] shows that 
DFT+U and GWU calculations on rutile using Uopt lead to very similar band gaps, 
2.83 eV and 2.85 eV, respectively. This suggests that the optimal Hubbard parameter 
f/opt determined for anatase is reasonably transferable to rutile Ti02. Similarly to 
the case of anatase Ti02, in the rutile polymorph the corrections to the individual 
eigenvalues calculated in DFT+f/ are ~0.5 eV higher than in GWU. 

In the case of anatase Ti02 we are not aware of published inverse photoemission 
data on single crystals of anatase, therefore it is not possible yet to perform a 
quantitative comparison between our calculated quasiparticle energies and experiment. 
We note however that the measured optical gap of anatase (3.2 eV [37]) is compatible 
both with our GW band gap (3.75 eV) and our GWU band gap (3.27 eV). In the case 
of rutile Ti02 the quasiparticle band gap has been measured in the range 3.3 ± 0.5 eV 
[38] . This indicates that our GW result (3.40 eV) and our GWU calculation starting 
from DFT+^7 (2.85 eV) are both compatible with experiment. The measured optical 
gap of rutile (3.0 eV [37]) is slightly larger than our calculated quasiparticle gap. 

A key result of the present work is that, in the case of anatase and rutile Ti02, 
the use of DFT+[/ as a starting Hamiltonian for calculating GW quasiparticle energies 
leads to band gaps which are 0.4 eV smaller than standard GW calculations starting 
from DFT ( [2^ [26] and table [1]). This relatively large difference between GW and 
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Table 1. Energies of critical points of anatase and rutile Ti02 calculated using various 
methods: plain density- functional calculations (DFT), calculations including Hubbard- 
like corrections with C/=7.5 eV (DFT-I-J7), and GqWq calculations starting from the 
DFT+f/ eigenstates and eigenvalues {GWU). For the sake of comparison with other 
studies [Ml 123 123 HH] we also report our results for GqWq calculations starting from 
DFT eigenstates and eigenvalues {GW). All values are in eV. The atomic positions 
and lattice parameters are taken from experiment |36j . In the case of anatase we also 
give the values calculated using our optimized lattice parameters between brackets. 
The choice of theoretical vs. experimental parameters lead to differences smaller than 
0.1 eV. 



Anatase TiO 
DFT 


2 


GW 




DFT+[/ 


GWU 




X24„ 0.00 
r25c 2.08 

Eg 2.08 


(0.00) 
(2.06) 
(2.06) 


-0.31 
3.39 
3.70 


(-0.30) 
(3.38) 
(3.68) 


0.48 (0.49) 
3.75 (3.70) 
3.27 (3.21) 


-0.14 
3.13 
3.27 


(-0.13) 
(3.07) 
(3.20) 


Rutile Ti02 
DFT 




GW 




DFT+[/ 


GWU 




r24. 0.00 

r25c 1.82 

Eg 1.82 




-0.19 
3.21 
3.40 




0.59 
3.42 
2.83 


0.08 
2.93 
2.85 





GWU raises the question of which Hamiltonian is a better starting point for many- 
body perturbation theory in strongly correlated systems. 

Before concluding we note that, since in the approach described here the optimal 
Hubbard parameter is determined through an iterative procedure, this scheme carries 
some similarities with self-consistent GW methods [8]. In the present work the self- 
consistent potential is constrained to be the DFT exchange and correlation potential 
plus an additional Hubbard-like term. It would be interesting to apply the same concept 
to the case of GW calculations starting from hybrid functional calculations [El 139] , 
where instead of the Hubbard U we would optimize the fraction a of Hartree-Fock 
exchange. 

In conclusion, we investigate the use of DFT+f/ as a starting point for performing 
many-body perturbation theory calculations on Ti02. We determine the Hubbard 
parameter of the non-interacting Hamiltonian in such a way that a GW calculation 
starting from the DFT+[/ band structure leaves the fundamental band gap unchanged. 
We find that the optimal Hubbard parameter of 7.5 eV yields the same band gap of 
3.27 eV in DFT-|-[/ and GWU . This value is considerably smaller than the value of 
3.70 eV obtained when applying GW corrections to the DFT band structure. The 
Hubbard parameter determined for anatase is transferable to the rutile polymorph 
of Ti02, and also in this case the calculated band gap is smaller than in standard 
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GW calculations starting from DFT. The difference between GW and GWU band 
gaps reported here calls for detailed photocmission experiments to establish whether 
DFT+^ represents a better starting point with respect to DFT for computing the GW 
quasiparticle energies of Ti02. 
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